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Abstract 

Background: Haemonchosis is a parasitic disease that causes severe economic losses in sheep industry. In recent years, 
the increasing resistance of the parasite to anthelmintics has raised the need for alternative control strategies. Genetic 
selection is a promising alternative but its efficacy depends on the availability of genetic variation and on the occurrence 
of favourable genetic correlations between the traits included in the breeding goal. The objective of this study was 
twofold. First, to estimate both the heritability of and the genetic correlations between growth traits and parasite 
resistance traits, using bivariate linear mixed animal models, from the phenotypes and genotypes of 1004 backcross 
lambs (considered as a single population), which underwent two subsequent experimental infestations protocols with 
Haemonchus contortus. Second, to compare the precision of the estimates when using two different relationship 
matrices: including pedigree information only or including also SNP (single nucleotide polymorphism) information. 

Results: Heritabilities were low for average daily gain before infestation (0.10 to 0.15) and average daily gain during the 
first infestation (0.1 1 to 0.16), moderate for faecal egg counts during the first infestation (0.21 to 0.38) and faecal egg 
counts during the second infestation (0.48 to 0.55). Genetic correlations between both growth traits and faecal egg 
count during the naive infestation were equal to zero but the genetic correlation between faecal egg count during the 
second infestation and growth was positive in a Haemonchus contortus free environment and negative in a 
contaminated environment. The standard errors of the estimates obtained by including SNP information were smaller 
than those obtained by including pedigree information only. 

Conclusions: The genetic parameters estimates suggest that growth performance can be selected for independently of 
selection on resistance to naive infestation. Selection for increased growth in a non-contaminated environment could 
lead to more susceptible animals with long-term exposure to the infestation but it could be possible to select for 
increased growth in a contaminated environment while also increasing resistance to the long-term exposure to the 
parasite. The use of molecular information increases the precision of the estimates. 



Background 

Haemonchus contortus {H. contortus) is a nematode that 
feeds on blood through the abomasal mucosa of bovine, 
ovine and caprine species [1]. The cost of H. contortus 
infection or haemonchosis for the production sectors of 
sheep farming in terms of anthelmintic treatments that 
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are currently the most popular control strategy, and the 
resulting economic loss have been estimated in different 
countries to be in the order of several million dollars per 
year [2-7]. Furthermore, anthelmintics tend to select the 
parasite population under treatment for resistance to the 
anthelmintic itself [8,9], which increases the cost of hae- 
monchosis even more. 

Both the long-term loss of efficacy and the growing 
public concern for the use of chemicals in food production 
fostered the research on alternative control strategies or 
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combinations of them [10], among which genetic selection 
is one of the most promising approaches [11-25]. Simula- 
tion studies based on evolutionary genetics [26] predict a 
breeding plans long-term outcomes and also the efficacy 
of genetic selection as a control strategy. Since genetic im- 
provement depends on the genetic parameters of the traits 
under selection, the estimates of these parameters must be 
as precise as possible for reliable long-term predictions. 
However, consistent estimates of the genetic correlation 
between production traits and parasite resistance traits 
have not been reported in the literature [27-33]. Since most 
of the estimates found in the literature are computed from 
observations in natural conditions, where it is not possible 
to precisely define neither the nutritional level of the diet 
nor the larval challenge on the pasture, the reason why no 
consistent estimates are available may be due to the inter- 
action between these two factors [34]. 

The first objective of this study was to estimate the 
heritability of average daily gain and faecal egg count 
from experimental observations, together with the gen- 
etic correlations between them. We report the results of 
an analysis performed on 1004 phenotypic records of 
growth traits and faecal egg counts collected on geno- 
typed (50 k SNPchip) back-cross lambs (25% Martinique 
black belly and 75% Romane), following two experi- 
mental infestations with H. contortus. The genetic pa- 
rameters have been estimated both by using pedigree 
information only and pedigree and SNP (single nucleo- 
tide polymorphism) chip information jointly. Computing 
the relatedness between individuals using pedigree infor- 
mation only is based on expectation and results in an es- 
timate corresponding to the average number of alleles 
shared by two individuals, for example: all the individuals 
belonging to the same full-sibs group would have a coef- 
ficient of 0.5 between each other, which means that it 
does not take into account the deviation from this aver- 
age caused by segregation and recombination. However, 
including molecular information makes it possible to 
compute the relatedness between individuals by identify- 
ing on a relatively dense map the actual number of alleles 
they share, which provides a more detailed estimate of 
the relatedness between individuals [35]. Since the obser- 
vations used in this study were collected on four large 
groups of half-sibs, the second objective of the study was 
to test whether including SNP information could help 
reach more precise estimates than using pedigree infor- 
mation only when the pedigree of the population is 
poorly informative. Although previous studies have already 
explored the amount of genetic variability for parasite 
resistance traits, this study features several novelties: the 
genetic parameters reported here are estimated from phe- 
notypes collected in experimental settings rather than nat- 
ural infestation, the growth traits analysed are the average 
daily gains before infestation and during infestation rather 



than the body weights, and finally the estimates reported 
here also feature molecular information rather than pedi- 
gree information only. 

Methods 

Experimental design 

The population in which the observations were collected 
from resulted from a back-cross mating scheme be- 
tween two pure-bred populations: Martinique Black 
Belly (MBB) and Romane (ROM). MBB is a tropical sheep 
breed, which is characterized by adaptation to heat-stress, 
to parasitism and to extensive raising conditions. The 
ROM sheep breed features good productive performances 
(both for meat production and prolificacy) and no selection 
for resistance to parasites. The pedigree used in the ana- 
lysis (Figure 1) was three generations deep and counted a 
total of 3164 animals. Four Fl sires were produced by 
crossing MBB and ROM individuals. The sires were mated 
by intra-uterine artificial insemination to 829 pure-bred 
ROM dams in order to obtain 1265 back-cross offspring 
(BC), the number of animals used in this study from each 
group of half-sibs were 282, 251, 247 and 223, respectively. 
All animals were managed as a commercial flock, the 
experimental protocol complied to the European ethical 
policy and was approved by the ethical committee "comite 
d'ethique CIRAD-INRA". 



MARTINIQUE 
BLACK 
BELLY 



ROMANE 




Figure 1 Scliematic representation of tlie pedigree structure. 

MBB: Martinique Blacl< Belly; ROM: Romane; Fl are the four Fl sires 
resulting from the cross between MBB sires and ROM dams; BC is the 
back cross population obtained by mating the Fl sires to ROM dams. 
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Genotypes 

A total of 1044 animals among the population of the 
back-cross lambs and their four Fl sires were geno- 
typed with the OvineSNPSO Beadchip (Illumina Inc., 
San Diego, CA). Quality control of the SNPs included 
the following tests: (1) 50 animals were genotyped 
twice in order to assess the technical reliability of the 
genotyping, which resulted in a value of 99.9%, (2) indi- 
viduals with a call rate below 98% and SNPs with a call 
rate below 97% or with a minor allele frequency below 
1% or featuring a deviation from expected heterozy- 
gosity or showing mendelian inconsistencies were dis- 
carded (p < 10' ) and (3) quality control of the genotypes 
resulted in 42 469 SNPs that comply with all quality 
checks. More details about the genomic information can 
be found in Salle et al. [36] . 

Phenotypes 

Phenotypes on growth traits and parasite resistance 
traits were collected on the BC animals only by perform- 
ing the following experimental protocol. The lambs were 
weaned around 64 or 45 days, depending whether their 
mother was either primiparous or not, and grew in a 
H. contortus free environment until the first experimen- 
tal infestation was performed. During this period the an- 
imals were kept in an H. contortus free environment and 
were weighed twice: at weaning and at the end of the 
growing period. The growing period ended with the 
beginning of the following experimental infestation 
protocol, which also determined when the environment 
was to be contaminated with H. contortus larvae: at 
around 90 days of age, the lambs received an oral inocu- 
lation of 10 000 L3 larvae of H. contortus (ENVT strain 
[37]) and around 41 days after the infestation, they re- 
ceived an anthelmintic treatment (LEVAMISOLE 5%, 
Vibrac S.A., Carros, France, 7.5 mg/kg live weight). Dur- 
ing the infestation, two faecal samples were collected, at 
25 and 35 days after infestation, and the animals were 
weighed on the day of treatment. Then, they entered a 
recovery period of 8 days, at the end of which they were 
infested again with the same infestation protocol. During 
the second infestation, two faecal samples were collected 
as before but animals were not weighed. During the 
whole protocol, the animals were fed ad libitum on a 
diet that covered largely their requirements. The faecal 
egg count in each sample was measured by a modified 
McMaster procedure [38]. The average of the two faecal 
egg count observations was computed for each infest- 
ation. The latter values were transformed by taking their 
fourth root in order to bring their distribution closer 
to normality. A further transformation was applied in 
order to scale the standard deviation to 1 and avoid 
zero values. The variables obtained were called: FECI 
(transformed faecal egg counts during the first infestation) 



and FEC2 (transformed faecal egg counts during the 
second infestation). The average daily gain from wean- 
ing to infestation (ADGO) and average daily gain during 
the first infestation (ADGl) were computed as follows: 



•^^""'^^ ■^"•^ ^ '^'^ ' — , where growthO 



ADGO = ^ and ADGl 



and timeO are the weight gain and the days running 
from the weaning day to the day of beginning of the 
first infestation, respectively; growthl and timel are the 
weight gain and the days running from the first day of 
the first infestation to the day of treatment of the first 
infestation, respectively. 

A transformation for scaling the standard deviation to 
1 and avoiding 0 s was applied on growth traits as well. 
The observations outside a range of 2.96 standard devia- 
tions around the average of each trait were considered 
atypical and excluded from the analysis. Finally, only the 
animals featuring a valid observation both on genotyping 
and on at least one trait were included in the analysis, 
which resulted in 40 animals being discarded and 1004 
animals being included. 

Statistical analysis 

Estimation of the genetic parameters was performed by 
considering the back-cross population as a single breed 
population. The breed proportions are taken into ac- 
count by the genetic relationship matrix in the model in- 
cluding SNP information but could not be taken into 
account in the model including pedigree information 
only due to convergence failure. The heritability of each 
trait and both the genetic correlation and phenotypic 
correlation between each pair of traits were estimated by 
bivariate animal mixed models, which were solved by 
the AIREML procedure implemented in AIRemlF90 soft- 
ware [39]. This software features by default the correction 
for the change in the definition of the base population so 
that the estimates obtained when using pedigree informa- 
tion only were comparable to those obtained when includ- 
ing molecular information [40,41]. The bivariate mixed 
model reads as follows: 



Xi 0 

0 X2 



bi 
b2 



Zi 0 

0 Z2 



where yi and y2 are the vectors of observations of trait 
one and two, respectively, Xi and X2 are incidence 
matrices relating each observation to its respective set of 
fixed effects and bi and \>2 are the vectors of the fixed 
effects: weight at weaning (for ADGO only) or weight at 
first infestation (for all the other traits), contemporary 
group (identified by year, season, weighting lot and 
infestation lot), sex and feeding mode (breast feeding 
or bottle feeding), ai and 82 are the vectors of random 
animals breeding values, with the associated incidence 
matrices Zi and Z2. Ci and 62 are the vectors of random 
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residuals. It is assumed that the random effects are nor- 
mally distributed and feature the following variance- 
covariance structure: 



V 
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where ffgu, 0^22 ^■^'^ ^i^i ^'^^ genetic variances 
and the genetic covariance between traits 1 and 2, 
CTg^j, 0^22 3nd (7e2j are the residuals variances and the 
residuals covariance between traits 1 and 2, 1 is an identity 
matrix and T is the genetic relationship matrix between 
the animals. 

The genetic and phenotypic correlations between each 
couple of traits were computed using two different gen- 
etic relationship matrices: one computed by using pedi- 
gree information only (PED model) and one computed 
by including both SNP and pedigree information jointly 
(SNPED model). The T matrix used in the PED model was 
computed according to Quaas [42] and without using mo- 
lecular information. The T matrix; used in the SNPED 
model was computed according to VanRaden [43], using a 
three-generation deep pedigree and the following weights: 
T = 0.95G + O.O5A22, where A22 is the relationship matrix 
between the genotyped animals computed by using the 
pedigree information only [42] and G is the genomic rela- 
tionship matrix among genotyped animals. In the software 
package used for this study [39], the G matrix is computed 
by default as follows: G = WK" ^W, [43] where W is a rect- 
angular matrix (number of animals by number of SNPs' 
alleles) with elements: w,y - 2py, where fij is a scalar 
equal to the number of copies of one allele an animal i has 
at locus j, Pj is the frequency of allele j in the population, K 
is the diagonal matrix of the scaling parameters with ele- 
ments: K^y = 2Z/?y(l -/?/)■ The weights of G and are used 
for bending the genetic relationship matrix and make it 
positive definite, as required for its inversion [43], this is 
similar to the bending procedure occurring in AIREML 
algorithms for keeping the variance covariance matrix 
positive definite [42] . 

The precision of the heritability and correlation esti- 
mates was computed by estimating their standard error 
according to the following formulas [44]: 




2C al,,a' 



^gii^pii 



where SE(hf) is the standard error of the estimate of the 
heritability of trait i, tr^,, and crj„ are the estimates of its 

/ gu pu f 0 \ 

genetic and phenotypic variances, respectively, Via^A, 
and C^ff^;;, (T^;;^ arc the variances of the 



V 
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estimated values and the covariance between the esti- 
mated values, respectively, obtained from the informa- 
tion matrix [44]: 



SE(rg) 



+ 



^c(a'^n,al22) c(a'^jj,ag2ij c{al22,'^g2i 



^gll^g22 



ffgn<^g2i 



(^g22'^g21 



where SE(r^ is the standard error of the estimate of the 
genetic correlation, rg is the estimated value of the gen- 
etic correlation; are ff^jj, (7^22 ^g2i the estimates of 
the genetic variance components described above and V(.) 
and C(..) are the variance of the estimates between paren- 
thesis and the covariance between the estimates between 
parenthesis, respectively. The same formula was used to 
compute the standard error of the phenotypic correlation, 
but by filling in the entries concerning phenotypic vari- 
ances and covariances. 

Significance tests 

The parameters under study are the ratio of two nor- 
mally distributed variables (heritability) and the ratio 
between a gaussian variable over the square root of the 
product of two gaussian variables (correlations). The 
sampling distribution of heritability can be approximated 
to a gaussian distribution under certain specific condi- 
tions only [45]. When these conditions are filled, the sig- 
nificance test for gaussian variables can be applied. 
However, the significance tests for the correlation coeffi- 
cient can be developed by deriving its confidence interval 
according to Fisher's Z-transformation [46]. Otherwise, 
both parameters can be tested by using a re-sampling pro- 
cedure such as delete-d jackknife: where d is the number 
of observations randomly discarded from the dataset and 
n^'^ < d < n (n is the total number of observations in the 
dataset) [47]. One thousand sub-samples of the whole 
dataset were created by randomly discarding 20% of the 
observations. Each parameter computed above was re- 
estimated from each sub-sample in order to build its 
empirical distribution. The empirical distribution of each 
parameter was used to compute the confidence interval of 
each parameter by taking its 2.5% quantile as the lower 
bound and its 97.5% quantile as the upper bound of each 
estimate. The null hypothesis "the estimate is not different 
from 0" was tested as follows: if the confidence interval of 
the estimate included 0, then the null hypothesis was not 
rejected, otherwise the alternative hypothesis "the estimate 
is different from 0" was accepted. In order to test whether 
the SNPED and SNP models converged to the same esti- 
mate, the distribution of the difference between the PED 
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and SNPED estimates was built for each estimate as 
follows: 

[di] = [Pedi-Snpedi] , 

where Pedi is the vector containing the n realizations of 
the estimate obtained from the PED model, Snpedi is 
the vector containing the n realizations of the estimate 
obtained from the SNPED model, di the vector of the 
differences between and each element of Pedi and 
Snpedi. The confidence interval of the distribution of 
the difference was computed as above. The null hypoth- 
esis "the difference between the estimate obtained from 
the PED model and the estimate obtained from the 
SNPED model is 0" was tested against the alternative hy- 
pothesis "the difference between the estimate obtained 
from the PED model and the estimate obtained from the 
SNPED model is not 0" as above as well. 

Results and discussion 

Phenotypic variation 

Table 1 shows the summary statistics of the variables 
analysed. The transformations applied to the raw faecal 
egg counts resulted in the profile of their distribution 
being closer to normality. The skewness and normalized 
kurtosis of FECI changed from 2.11 to -0.44 and from 
6.56 to 0.53, respectively; the skewness and normalized 
kurtosis of FEC2 changed from 3.87 to 0.14 and from 
23.92 to -0.74, respectively. The number of observations 
on each trait together with the average, standard devi- 
ation, minimum and maximum of the raw observations 
are in Table 1. ADGl was significantly lower than ADGO 
(p_value < 0.0001), indicating that infested animals had a 
slower growth than the parasite-free animals, as ex- 
pected due to the infestation [1]. 

Table 2 shows the estimate of the phenotypic correla- 
tions (below the diagonal) obtained from the two models 
for each pair of traits. Although the SNPED and PED 
models did not always converge on the same value, ac- 
cording to the significance tests described above, these 
estimates were not significantly different between the 
two models. The estimates of the phenotypic correla- 
tions between ADGO and both FECI and FEC2 were not 
significantly different from 0: -0.01 (SE = 0.15) and 0.01 



(SE = 0.18) for the PED model and 0.02 (SE = 0.11) and 
0.04 (SE = 0.11) for the SNPED model. These results 
suggest that the phenotype for growth rate in a H. con- 
tortus free environment was unrelated to the parasite 
resistance phenotype. However, the estimates of both 
the phenotypic correlations between ADGl and FECI 
and between ADGl and FEC2 were negative: -0.24 
(SE = 0.15) and -0.20 (SE = 0.19) for the PED model 
and -0.23 (SE = 0.11) and -0.19 (SE = 0.11) for the SNPED 
model. These results suggest an inverse proportionality 
between the growth rate and the parasite burden, in ac- 
cordance with the finding that contaminated animals had 
a slower growth than non-contaminated animals. The 
average faecal egg count during the second infestation was 
significantiy lower than the faecal egg count during the 
first infestation (p value < 0.0001), which suggests that the 
development of a specific immune response was triggered 
by the first infestation, that enhanced the intrinsic resist- 
ance of the animals to subsequent infestations [48]. 
Furthermore, the positive estimate of the phenotypic cor- 
relation between FECI and FEC2, 0.46 (SE = 0.43) for the 
PED model and 0.62 (SE = 0.20) for the SNPED model, 
shows that the animals featuring higher (or lower) than 
average FECI are likely to express higher (or lower) than 
average FEC2, and vice versa. This suggests that a re- 
peatable variation in susceptibility occurs within the 
population. 

Genetic variation 

Table 2 summarizes the estimates of the heritabilities of 
each trait (block diagonal) and the estimates of the gen- 
etic correlations (above the diagonal) between each pair 
of traits obtained with the PED and SNPED models, to- 
gether with the standard error of each estimate (between 
brackets). Due to the pedigree structure that includes 
only four sires, both the standard errors and the 95% 
confidence intervals of the estimates were indeed large 
(in particular, those of the PED model), which led to no 
significant difference between the estimates obtained 
from the two models. The results obtained were in gen- 
eral coherent between models, except for the genetic 
correlation between ADGO and FECI, which was posi- 
tive with the SNPED model and negative with the PED 
model. However, the latter estimate had a standard error 



Table 1 Descriptive statistics of the raw observations 



Raw observations 


Number of observations 


Average 


Standard deviation 


Min 


Max 


ADGOraw (g/day) 


997 


293.3 


694 


93.0 


504.5 


ADGl RAW (g/day) 


963 


102.8 


43.3 


-27.0 


230.0 


FECI RAW (eggs/g) 


987 


10494 


9827 


0 


75898 


FEC2raw (eggs/g) 


967 


2724 


4259 


0 


42667 



ADGOraw = average daily gain before infestation; ADGl raw = average daily gain during tlie first infestation; FECI raw = faecal egg count during the first infestation; 
FEC2raw = faecal egg count during the second infestation. 
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Table 2 Heritabilities, phenotypic correlations and genetic 
correlations obtained by the SNPED model and the PED 
model 



Traits 


Model 


ADGO 


ADGl 


FECI 


FEC2 


ADGO 


SNPED 


0.15 (0.07)* 




0.11 (0.47) 


0.57 (0.38)* 




PED 


0.10(0.08)* 




-0.52 (1.06) 


0.25 (0.85) 


ADGl 


SNPED 




0.11 (0.06)* 


-0.12 (0.58) 


-0.54 (0.53)* 




PED 




0.16 (0.04)* 


-0.19 (0.80) 


-0.48 (0.67)* 


FECI 


SNPED 


0.02 (0.11) 


-0.23 (0.11)* 


0.38 (0.04)* 


0.62 (0.20)* 




PED 


-0.01 (0.15) 


-0.24(0.15)* 


0.21 (0.06)* 


0.46 (0.43)* 


FEC2 


SNPED 


0.04 (0.11) 


-0.19(0.11)* 


0.31 (0.08)* 


0.48 (0.06)* 




PED 


0.01 (0.18) 


-0.20 (0.19)* 


0.29 (0.14)* 


0.55 (0.09)* 



ADGO = average daily gain before infestation; ADG1 = average daily gain 
during the first infestation; FECI = faecal egg count during the first infestation; 
FEC2 = faecal egg count during the second infestation; SNPED refers to estimates 
obtained by using the joint pedigree and molecular information relationship 
matrix (SNPED model); PED refers to the pedigree-only relationship matrix (PED 
model); the correlations between ADGO and ADGl were much more sensitive 
than others to the starting values used for the estimation and to re-sampling, and 
are not presented; heritabilities are on the block diagonal, genetic correlations 
are above the diagonal and phenotypic correlations are below the diagonal; 
standard errors of the estimates are between brackets; superscript * marks the 
estimates that were significantly different from 0. 

as large as half the parameter space, which resulted in 
the confidence of the difference between the two esti- 
mates to include 0. The reason for this possible incon- 
sistency cannot be defined by the data available for this 
study. We can only speculate that it could be the result 
of the segregation variance captured by the SNP chip, 
because the phenotypes analysed were collected on the 
back-cross offspring of only four sires. This pedigree 
structure is indeed poorly informative if the estimate is 
computed by using pedigree information only, which is 
clearly shown by the huge standard error obtained with 
the PED model. However, this does not explain why 
such a big difference in the estimates occurs between 
some pairs of traits only. The heritability of ADGO was 
low for both models, 0.10 (SE = 0.08) for the PED model 
and 0.15 (SE = 0.07) for the SNPED model, and is 
close to the value 0.17 found in the review of Safari 
et al. [27]. The heritability of ADGl was low as well, 0.16 
(SE = 0.04) and 0.11 (SE = 0.06) for the PED and the 
SNPED model respectively and no estimates were found 
in the literature for ADGl. The estimates for the herita- 
bilities of both faecal egg counts were found higher than 
those obtained in previous studies [29-32]: FECI was 
0.21 (SE = 0.06) and 0.38 (SE = 0.04) for the PED and the 
SNPED model, respectively; FEC2 was 0.55 (SE = 0.09) 
and 0.48 (SE = 0.06) for the PED and the SNPED model, 
respectively. The reason for finding a higher value than 
in previous studies could be twofold. First, it could be 
due to the fact that most of the estimates found in the lit- 
erature are computed from observations in natural rather 
than experimental conditions, the latter of which allows 



controlling more strictly the environmental conditions 
and hence could reduce phenotypic variation. The sec- 
ond reason could be the occurrence of breed specific 
alleles that segregate within the back-cross population, 
which inflate the genetic variance compared to a pure 
breed population. These values confirm the availability of 
a moderate genetic potential in sheep that could be 
exploited to enhance resistance to parasites. 

Concerning previous estimates of the genetic correl- 
ation between growth traits and faecal egg counts, no 
other estimations of these parameters based on experi- 
mental infestation were found in the literature. Further- 
more, previous studies on similar traits based on natural 
infestations do not show consistency among them [29-32], 
which could be explained by the uncontrolled variation in 
the larval challenge, in the pathogenicity of the parasite in 
each population, in the feed intake and the interaction 
between these three factors [34]. 

Concerning the genetic correlations between growth 
traits (ADGO and ADGl) and parasite resistance during 
the naive infestation (FECI), the following picture can be 
drawn. In accordance with the estimates of the pheno- 
typic correlation between ADGO and FECI, the genetic 
correlations between these traits were also not signifi- 
cantly different from 0: -0.52 (SE = 1.06) for the PED 
model and 0.11 (SE = 0.47) for the SNPED model. The 
same results were obtained with the estimate of the 
genetic correlations between ADGl and FECI (although 
their phenotypic correlations were negative according 
to both models): -0.19 (SE = 0.80) for the PED model 
and -0.12 (SE = 0.58) for the SNPED model. These re- 
sults suggest that during the naive infestation, the geno- 
type for growth (ADGO and ADGl) could be expressed 
independently from the genotype for parasite resistance 
(FECI). 

The results obtained for the genetic correlations be- 
tween growth traits and the long-term resistance (FEC2) 
were on the contrary significantly different from 0. On 
the one hand, the genetic correlation between ADGO 
and FEC2 was positive but not significantly different 
from 0 according to the PED model (0.25, SE = 0.85) and 
positive (0.57, SE = 0.38) according to the SNPED model. 
According to the approximate standard errors, the posi- 
tive estimate obtained with the SNPED model is more 
reliable and suggests that if animals were selected for 
growth in a parasite-free environment, a correlated selec- 
tion response for lower long-term resistance to gastro- 
intestinal parasites could occur as well. This estimate 
supports the hypothesis that enhancing growth traits 
could come to a cost to the sheep's long-term suscep- 
tibility to parasite infestations, and vice versa [49,50]. 
On the other hand, the correlation between ADGl and 
FEC2 was consistently negative between models: -0.48 
(SE = 0.67) for the PED model and -0.54 (SE = 0.53) for 
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the SNPED model, which suggests that growth during the 
naive infestation can be enhanced together with long-term 
resistance to the infestation within a single purebred line. 

The genetic correlation between FECI and FEC2 was 
0.46 (SE = 0.43) and 0.62 (SE = 0.20) for the PED and 
SNPED models, respectively, which suggests that these 
traits have different determinisms. While FECI repre- 
sents a measure of the parasite resistance expressed by a 
naive lamb, FEC2 is a measure of the parasite resistance 
expressed by an immunized lamb, and indeed the mech- 
anisms by which these types of animals respond to the 
infestation are different [48]. 

The estimate of genetic correlation between ADGO 
and ADGl was unstable due to its sensitivity to the starting 
values used for its estimation and is not reported. 

Standard errors 

Table 3 shows both the ratio of the standard errors of 
the estimates obtained from the PED model over the 
standards error obtained from the SNPED model, which 
ranged from 1.04 to 2.25 and also the ratio between the 
width of the confidence intervals of the estimates ob- 
tained from the PED model over the width of the confi- 
dence intervals obtained from the SNPED model, which 
spanned an interval between 0.93 and 4.21. According to 
the ratio of the standard errors, the SNPED model al- 
ways converged to more precise values, while, according 
to the ratio of the width of the confidence intervals, the 
estimate of the phenotypic correlation between growth 
traits and FECI obtained from the PED model was 
slightly more precise. 



Table 3 Ratios between the precision estimators for the 
PED model and for the SNPED model 



Traits 


ADGO 


ADGl 


FECI 


FEC2 




^^1.12 






^^2.25 


ADGO 


°142 




•='1.26 








"'4.21 






=^^1.23 




=^1.26 


ADGl 














^'2.37 


°1.30 


"1.57 




^^1.37 


^^1.33 


=^^172 


^^2.22 


FECI 


"0.93 


^'0.98 


°1.5 


"'2.2 




=^^1.73 


^^1.68 




=^1.67 


FEC2 


°141 


"'1.30 


"='1.25 






"'1.9 



ADGO = average daily gain before infestation; ADG1 = average daily gain 
during the first infestation; FECI = faecal egg count during the first infestation; 
FEC2 = faecal egg count during the second infestation; the table shows both 
values of the ratio between the approximate standard errors (marked with 
superscript SE) and the ratios of the width of the confidence intervals 
obtained from the empirical distribution {marked with superscript CI) of the 
parameters estimates obtained by using either pedigree relationship matrix or 
the joint pedigree and molecular information relationship matrix; the ratios on 
the diagonal refer to the heritability estimates; the ratios above the diagonal 
refer to the genetic correlation estimates; the ratios below the diagonal refer 
to the phenotypic correlation estimates; the correlations between ADGO and 
ADG1 were much more sensitive than others to the starting values used for 
the estimation and to re-sampling, and are not presented. 



The results obtained show that including SNP infor- 
mation in the computation of the relationship matrix 
between individuals can increase the precision of the 
genetic parameter estimates up to twice the precision 
obtained by using pedigree information only [51]. The 
increase in precision can be explained by the fact that 
SNP information allows to compute more precisely than 
pedigree information what proportion of genome two in- 
dividuals actually have in common. The pedigree struc- 
ture in the data available for this study was not ideal to 
estimate genetic parameters by pedigree information 
only because all observations are recorded on a popula- 
tion of animals composed of four groups of half-sibs. 
Such a structure causes the pedigree-based relationship 
matrix to predict that within each group of half-sibs all 
animals share one quarter of the sire's genome. Whereas, 
the marker-based relationship matrix allows capturing 
the segregation variance, which means capturing the 
random deviation of the proportion of genes shared by 
two individuals around the expected proportion of shared 
genes according to the pedigree [35]. 

Conclusions 

According to the results obtained by the model includ- 
ing both pedigree and molecular information, the geno- 
types for growth and for resistance to naive infestation 
can be selected for independently. However, the genetic 
correlations between long-term parasite resistance traits 
and growth traits were different from 0 and suggest that 
increasing growth performance in a H. contortus free 
environment could result in more susceptible animals, 
whereas growth performance in a contaminated environ- 
ment can be increased while enhancing long-term resist- 
ance to H. contortus. The two results taken together can 
also be interpreted as an indication of genotype by envir- 
onment interaction affecting growth expressed across the 
two environments [52]. The model that includes pedigree 
information only converged to similar results, except for 
the genetic correlation between growth before infestation 
and faecal egg count during the first infestation which was 
affected by a very large standard error. The reason for this 
inconsistency needs further investigation. 

This study shows that, when the pedigree is poorly 
informative using molecular information and pedigree 
information jointly result in more precise genetic param- 
eters than using pedigree only. 
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